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Single-layer and Bilayer of graphene are new classes of two-dimensional electron systems with 
unconventional band structures and valley degrees of freedom. The ground states and excitations in 
the integer and fractional quantum Hall regimes are investigated on torus and spherical geometries 
with the use of the density matrix renormalization group (DMRG) method. At nonzero Landau 
level indices, the ground states at effective filling factors 1, 1/3, 2/3 and 2/5 are valley polarized 
both in single-layer and bilayer graphenes. We examine the elementary charge excitations which 
could couple with the valley degrees of freedom (so called valley skyrmions). The excitation gaps 
are calculated and extrapolated to the thermodynamic limit. The largest excitation gap at effective 
filling 1/3 is obtained in bilayer graphene, which is a good candidate for experimental observation 
of fractional quantum Hall effect. 
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I. INTRODUCTION 



Hamiltonian for the K-vallcy of the SLG is written as 



A recent experimental realization of single-layer 
graphene (SLG) sheets^ has made it possible to confirm 
a number of theoretical predictions of intriguing electric 
properties of massless Dirac fermion systems,— including 
unconventional quantum Hall effects (QHE)^ii with the 
half-integer Hall conductivity^ 
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aX V — ±2, ±6, ±10, • • • , where a factor 4 is the Landau 
level (LL) degeneracy, resulting from spin and valley (re- 
ferred to K and K') symmetry in graphene. v — 2iTpgp 
is the filling factor, £b = \J h/eB is the magnetic length, 
p is the carrier density measured from the charge neutral 
Dirac point. 

As many phenomenological insights were found in 
the quantum Hall regime of SLG^>ii^i^ii^iaiiaiiiii^i^iiii^, 

bilayer graphene (BLG) potentially exhibits rich 
physicSi ^^'^^ Experiments showed that the Hall conduc- 
tivity of unbiased graphene bilayer in strong magnetic 
fields is given by 
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with |n| > 1.— The absence of (Jxy ~ plateau and the 
double height jump of the Hall conductivity between v = 
—4 and 4 indicate eight-fold degeneracy at the neutrality 
point. 

These quantization rules of the Hall conductivity for 
SLG and BLG originate with the characteristic energy 
spectrum in a magnetic field. The low-energy band struc- 
ture in SLG consists of Dirac cones located at the in- 
equivalent Brillouin zone corners K and K'. Using the 
magnetic ladder operator a = (tt^ — iTTy)lB / ^f^h, where 
TT = p — eA is the kinetic momentum, the single-particle 
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and Ti^^'^'' is the transpose of TC^^^^' . Where vp is the 
Fermi velocity of SLG, 2x2 matrices act in the sublat- 
tice degrees of freedom in graphene. The eigenenergy of 
eq.© is given by 
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The eigenvector is jO)^^'^'^'' 
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K = (0, |0))* for n = 0, and 
l),±|n))7^/2 for n > 1. Here \n) is 
fltfl with an 
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the eigenvector of the number operator N ; 
eigenvalue n. 

In contrast to SLG, BLG has an ordinary parabolic 
spectrum in the vicinity of the neutrality point. In a mag- 
netic field, the effective Hamiltonian for BLG is written 
in the forms^i 
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and its eigenenergy is 



En = ±huJc\Jn{n - 1). 
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Here lJc ~ eB/m with m being the effective mass of 
BLG. The eigenvector of eq.l^ is |0,0)^^°^ (0, |0))*, 
(0, |1))* for the zero energy level, and 
(|n-2),±|n))7A/2 forn> 2. 
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Recent experimental studies on SLG in a sufficiently 
strong magnetic field revealed new quantum Hall states 
at = 0, ±1, ±4^ where the electron-electron interaction 
may play a crucial role. Here relevant energy scales in 
graphene in a magnetic field are 

(i) LL separation around the neutrality point, 
V2hvp/iB ^ 400 V-B[T] [A^ for SLG, while V2hLJc ^ 
30 X (B[T])[K] for BLG. 
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gfiB\B\c,1.5x{B[T])[K], 



(ii) Zeeman coupling, 
and 

(iii) The Coulomb energy, e'^/ds ^ lOOy/ B[T] [K]. 

The activation energy measurements at these addi- 
tional QHE states in SLG^ have shown that si v = ±4 
the gap has linear B dependence and reasonably cor- 
responds to A^, indicating Zeeman spin splitting. At 
v = ±1, on the other hand, the gap is approximately 
scaled by a/B, that indicates the gap originates from 
the Coulomb interaction. The latter behavior is consis- 
tent with the quantum Hall ferromagnetism (QHF}i2, in 
which valley degrees of freedom, referred as pseudospins, 
spontaneously split via the exchange energy at all inte- 
ger fillings Based on the above scenario, it 
has been anticipated that interactions will drive quantum 
Hall effects also in BLG, at the octet's seven intermedi- 
ate integer filling factors when magnetic field is strong 
enough or disorder is weak enough.— 

In this work, we focus on the many-body states in 
nonzero LLs of SLG and BLG. We introduce the effective 
filling factors of the topmost partially filled nth LL, 
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for BLG. Here < 4, and n > 2 for BLG. 

The ground state at Vn = \ is fully spin and valley 
polarized, and the wave function can be represented by 
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where we assign the valley K and K' in graphene as 
z-component of pseudospin r = K or K'. We omit 
real spins, which are assumed to be fully polarized by 
the strong Zeeman splitting. Excitations from the sym- 
metry broken states are described by (pseudo)spin wave 
and (pseudo)spin textures called skyrmion o^^i^^'^^i^^'^"^ or 
other types, depending on the LL index n4 

In the previous work, we have studied the fractional 
quantum Hall ferromagnetic states in SLG at Un ~ 1 
and l/3.i^ As a consequence of the relativistic nature 
of electrons in SLG, the effective electron-electron in- 
teractions in n 7^ LLs differ from that in conven- 
tional two-dimensional systems, while the n = LL is 
equivalent liiiSiii We have neglected the spin degrees of 
freedom of electrons by assuming strong Zeeman cou- 
pling, although valley degrees of freedom have been taken 
into account since there is no external symmetry breaking 
field. We have shown that the ground states at i/„ = 1 
and 1/3 in n = and 1 LLs of SLG are fully valley 
polarized, while elementary charge excitations consist of 
pseudospin-singlet, namely valley-skyrmions ,20,21,22,23,24 

In this paper we extend the quantum Hall ferromag- 
netism to BLG and numerically show that the ground 
states at = 1 and 1/3 in the n ~ 2 LL of BLG are 



fully valley polarized. In contrast to the fact that ele- 
mentary charge excited states are pseudospin-singlet at 
1^2 = 1, those at = 1/3 of BLG are pseudospin polar- 
ized. Namely, the Laughlin type quasiparticle (quasihole) 
excitations dominate over skyrmion type. As a conse- 
quence, the charge gap at 1^2 — 1/3 in BLG is almost 
twice larger than that at I'd = 1/3 in SLG. Therefore 
clean BLG samples arc better candidates to observe frac- 
tional QHE than SLG. We also study = 2/3 and 2/5 
states in SLG and BLG. We claim that the ground at 
i^n = 2/3 and 2/5 states are pseudospin singlet in the 
n ~ LL of SLG, while fully pseudospin polarized in the 
n=lLL of SLG and the 71 = 2 LL of BLG. These resuffs 
are summarized in Table I. 



II. MODEL AND METHOD 

We start with the projected Coulomb interaction 
Hamiltonian onto a certain LL, and study charge and 
valley excitations, where we treat the valley degrees of 
freedom K and K' in the language of the pseudospin, 
while real spin degrees of freedom are supposed to be 
frozen by the Zeeman splitting. The LL mixing is also 
neglected in the following. We calculate the exact wave 
function of the ground state and low energy excited states 
in SLG and BLG, basing on the density matrix renormal- 
ization group (DMRG) method^i^^i^l and examine the 
existence of the skyrmion excitations in valley degrees of 
freedom not only at integer fillings = 1 but also at 
fractional fillings f„ = 1/3, 2/5 and 2/3 for the LL in- 
dices n = 0, 1 and 2. Charge gaps at these fractions are 
extrapolated to the thermodynamic limit. For our pur- 
pose the DMRG method is quite useful, since it needs to 
treat a large number of basis of the many-body Hilbert 
space when the pseudospin degrees of freedom are intro- 
duced in the fractional QHE systems. Taking account of 
valley degrees of freedom, we apply the DMRG method^^ 
on torus2Si21 and spherical geometries i^^i^^i^^i^°i'^^i^^ 
The projected Hamiltonian onto the nth LL is written 
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where is the guiding center coordinate of the ith par- 
ticle. The relativistic form factor in the nth LL of SLG 
is written as^ 



(11) 



and 



FLfiil) = (1/2) [L\ni (gV2) (9V2)] . (12) 

For the n(> 2)th LL in BLG, 

Filf\q) = (1/2) (gV2) + i|„|-2 (gV2)] . (13) 
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Here are the Laguerre polynomials. In the spher- 

ical geometry, it is convenient to write the Hamiltonian 
as 
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where Pij [m] projects onto states in which particles i 

and j have relative angular momentum hm, and T^m"'' 
is their interaction energy in the nth LLi^ Using above 
form factors, the pseudopotential o^^'^^ are given by 

Kl"^ = r ^qV{q)e-'^\F^{q)fL„M^). (15) 

The corresponding integrals for electrons on the surface 
of a sphere which are used in the present work are de- 
scribed in refs. I29ll30ll3l] . Note that our Hamiltonian 
eqs. (fTO|) and (|14p have SU(2) symmetry in the valley 
degrees of freedom. A symmetry breaking correction to 
eqs. pO| and p4|) which stems from the lattice structure 
of SLG^iiS and BLG is order of a/ls [a being a lattice 
spacing) in units of j els and neglected in the following. 

We calculate the ground state wave function using the 
DMRG method^^i^ which is a real space renormaliza- 
tion group method combined with the exact diagonaliza- 
tion method. The DMRG method provides low-energy 
eigenvalues and corresponding eigenvectors of the Hamil- 
tonian within a restricted number of basis states. The 
accuracy of the results is systematically controlled by 
the truncation error, which is smaller than 10~^ in the 
present calculation. We investigate systems of various 
sizes with up to 40 electrons in the unit cell keeping 1400 
basis in each blockj ^^'^^ 

The sphere geometry is useful to extrapolate energy 
gaps to the thermodynamic limit. In the sphere geom- 
etry, the pseudospin (valley) polarized ground state at 
Vn = 1/q {q being an odd integer), the Laughlin state?^ 
realizes when the total flux Ncj, is given hy^ 
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where is the number of electrons in the system. Ele- 
mentary charged excitations from this pseudospin polar- 
ized ground state correspond to the ground state config- 
urations of the system with additional/missing flux ±1. 
At Vn ~ 1/q, we study two types of excitations: Laugh- 
lin's quasiholcs (quasiparticles)^ and skyrmion quasi- 
holes (quasiparticles)i^i2ii22ii^i2i Laughlin's quasiholes 
(quasiparticles) correspond to pseudospin polarized ex- 
citations with ±1 flux, whose creation energy is given 

by 

Af =EiN^±l,P=l)-EiN^,P^l), (17) 

where ± represents quasiholes and quasiparticles, respec- 
tively, and P is the polarization ratio of the pseudospin, 
i.e. P = {Nk - Nk')/{Nk + Nk') with Nk {Nr') being 
the number of electrons in K {K') valley. 
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FIG. 1: The lowest charge excitation gaps in SLG and BLG 
at i/n = 1 in the n = 0, 1 and 2 LLs in the spherical ge- 
ometry. Closed circles represent the pseudospin unpolarized 
P = (skyrmion) excitation gaps and open circles represent 
nearly pseudospin-polarized P = 1 — 2/Ne excitation gaps. 
The crosses on the vertical axis represent results obtained by 
Hartree-Fock calculations.— 



Skyrmion quasiholcs (quasiparticles) correspond to 
pseudospin singlet excitations, and their creation energy 
is given by 



E{N^±1,P = Q) - E{N^,P ^1) 



(18) 



which could be smaller than 

The activation energy, referred as the gap in the fol- 
lowing, is given as a sum of these quasiholc and quasipar- 
ticle energies, Ac = A+ -I- A~ for pseudospin polarized 
(Laughlin-like) excitations, and A^ = A+ + A^ for pseu- 
dospin unpolarized (skyrmion-likc) excitations. 

At Un = 2/5 and 2/3, the pseudospin-polarized ground 
state and unpolarized ground state compete each other. 
Unfavorably, in the spherical geometry, these ground 
states realize at different configurations (the number of 
electrons and the total flux). Precisely, the pseudospin- 
polarized Vn = 2/5 state with electrons occurs when 
the total flux is given by = (5/2)iVe - 4. On the 
other hand, the pseudospin-singlet (unpolarized) state is 
realized at = (5/2)7Ve — 3. Since the finite size ef- 
fects are different between these different configurations, 
it is difficult to study the pseudospin-polarizability in the 
ground state in the sphere geometry. To avoid this diffi- 
culty, we utilize the torus geometry in which two ground 
states with P = and P = 1 are realized in the same 
configuration: = (5/2)7Ve. We use the sphere ge- 
ometry to extrapolate the energy gaps to the thermody- 
namic limit. The f„ = 2/3 states are studied as well in 
the following section. The elementally charge excitations 
at f„ — 2/3 are obtained by adding or removing single 
electron in the system with ±1 flux. 
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l/Ne 

FIG. 2: Expectation values, (cJ^rCmr)) in the lowest pseu- 

dospin (valley) unpolarized excited state with extra one flux FIG. 3: The pseudospin (valley) polarized excitation gap Ac 
at (a) = 1 in BLG and at (b) fi =2/5 in SLG. and the pseudospin (valley) unpolarized (skyrmion) excitation 

gap As at = 1/3 in SLG {n = and 1 LLs) and in BLG 

(n = 2 LL). 



III. RESULTS 

A. Vn = 1 and 1/3 states 

At t'n = 1 integer fillings, the pseudospins of valley 
degrees of freedom in SLG and BLG are completely po- 
larized by the exchange Coulomb interaction. This is es- 
sentially the same as the case of usual quantum Hall sys- 
tems with spin degrees of freedom. The fully pseudospin 
polarized P = 1 quasiparticle excitations from the above 
ground state need energy to the next higher LL, and such 
excitations have large gap compared with the unpolarized 
P = or partially polarized < P < 1 excitations in the 
same LL. We therefore calculate the unpolarized and par- 
tially polarized excitations by the DMRG method, and 
study the pseudospin polarization of the lowest charge 
excitation and the gap in the thermodynamic limit. 

The pseudospin polarization and the energy of the low- 
est charge excitation are shown in Fig. 1 for various sizes 
of system in spherical geometry. This figure shows that 
the pseudospin unpolarized (P = 0) state is the lowest 
excited state when the system size is large enough in SLG 
(n = 0, 1, and 2) and BLG (n = 2). The extrapolated 
value of the pseudospin unpolarized excitation gaps 
in units ofe'^/{e£B) are 0.63, 0.28 and 0.42 in the ti = 0, 1 
and 2 LLs of SLG, and 0.48 in the n = 2 LL of BLG. 

In the n = 2 LL, the lowest quasiparticle excitation in 
small systems (iVe ^ 20) has large pseudospin polariza- 
tion P = l-2/Ne both for SLG and BLG. These resuhs 
show instability of the unpolarized P = excitations in 
small systems. We find the first order transition in the 
lowest charge excited state when the number of electrons 
A^e exceeds 18 in SLG and 24 in BLG. This is consistent 
with the expectation that the unpolarized excitations are 
unstable in higher LLs because of the long-range nature 
of the effective exchange interaction, which increases en- 
ergy of skyrmion-like pseudospin unpolarized state. 



To study the pseudospin structure in the unpolarized 
excited states, we compare our numerical results with 
the Hartree-Fock (HF) trial states of skyrmions. The HF 
trial state of quasihole skyrmions at = 1 is written in 
the formii^ 

l*.fc>= n [^m(^lnK+Pmcl + ,KM, (19) 
jn=-N^/2 

where (c|,^^c„^) = and (cj„^,c„^,) = 

we have calculated the expectation values {c\^^c^^) from 
the wave function obtained in the present DMRG study. 
The results for U2 = 1 in BLG indicate they are approx- 
imately given by (cj„^c„^) = 1/2 =F m/N^ for t = AT 
and T = K' , respectively as shown in Fig. 2 (a). Similar 
results are also obtained for = 1 and 1 /3 in the n = 
and 1 LLs of SLG, that shows the pseudospin unpolar- 
ized elementally charge excitations in BLG at = 1 are 
pseudospin (valley) skyrmions. 

The elementally charge excitation gaps at fractional 
fillings i^n = 1/3 are shown in Fig. 3. The pseudospin 
(valley) polarized P = 1 excitation gap Ac for ti = 
SLG is 0.101 in units of e^/(e£B) in the thermodynamic 
limit, which is in a good agreement with the previous 
work,^ In the n = 1 LL of SLG, A^ is 0.115, which 
is larger than 0.101 in the n = LL of SLG. In BLG, 
Ac = 0.103 in the n = 2 LL, which is also slightly larger 
than that in the n = LL of SLG. This enhancement 
of the Ac in higher LLs is consistent with the increase 
of the difference of Haldane's pseudopotentials between 
m = 1 and 3; - Kj" is 0.168 in units of eV(e^B) for 
BLG in the 71 = 2 LL and 0.198 for SLG in the n = 1 
LL, while it is 0.166 for SLG in the n = LL. 

As seen in conventional 2DEGs, v — 1/3 (valley or 
spin) unpolarized excited states can be lower in energy 
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FIG. 4: Size dependence of the energy difference between 
the fully pseudospin polarized P = 1 state i?p=r and the 
unpolarized P = states Ep^o at = 2/5 in SLG in a torus 
geometry. Inset shows size dependence and V2 dependence of 
Ep=i ~ Ep=o at u„ = 2/5 in SLG with n — 1. 



FIG. 5: The pseudospin (valley) polarized excitation gap A^ 
and the pseudospin (valley) unpolarized excitation gap As at 
Vn — 2/5 in n — and 1 LLs of SLG and in the n = 2 LL of 
BLG. 



than the polarized excited states. Indeed, as shown in 
Fig. 3, the unpolarized excitation gap As in SLG is 0.05 
in the n = LL and 0.03 in the n = 1 LL, which are much 
smaller than the polarized excitation gap in the n ~ 
and 1 LLs, respectively. In the case of BLG, however, we 
find that valley unpolarized excited states have higher 
energy than polarized excited states. Consequently, the 
activation energy at = 1/3 in BLG is 2-3 times larger 
than that in SLG, that indicates an advantage for the 
observation of fractional QHE in BLG. 

B. — 2/5 and u„ — 2/3 states 

According to the composite fcrmion theory^, v = 
to/(1 + 2m) fractional quantum Hall states arc mapped 
onto the I'oft = m integer quantum Hall states. In con- 
trast to the case of 1/ = 1/3, where the ground state 
is mapped onto the n = 1 pseudospin polarized inte- 
ger quantum Hall state, 1^ — 2/5 fractional quantum 
Hall state is mapped onto the n = 2 integer quan- 
tum Hall state, where the lowest LL is doubly occupied 
by pseudospin-up and pseudospin-down electrons. We 
therefore expect pseudospin unpolarized ground state at 
f„ = 2/5. Here we calculate the ground state pseudospin 
polarization of SLG and BLG at = 2/5 to see whether 
this naive expectation is correct even in graphene. 

The pseudospin-polarization in the ground state is 
studied by calculating the polarization energy i?p=i — 
Ep^o which is the energy difference between the polar- 
ized P = 1 state and the unpolarized P = state. This 
energy difference is plotted in Fig. 4 for the n = and 1 
LLs in SLG as a function of TVe- This figure clearly shows 
that the ground state pseudospin at f„ = 2/5 in the 



n = LL of SLG is unpolarized, because i?p=i — £'p=o 
increases with the increase in the number of electrons 
Ne. This is consistent with the prediction of the com- 
posite fermion theory. On the other hand, in the n = 1 
LL of SLG, Ep=i — Ep=o decreases with the increase in 
Ne. This behavior suggests the polarized ground state. 

To confirm the polarized ground state at i'™ = 2/5 in 
the n = 1 LL of SLG, we slightly change the Haldane's 
pseudo potential V2, which acts only for electron pairs 
whose relative angular momentum m is 2. Since anti- 
symmetrizcd wave function of polarized pseudospin state 
does not contain electron pairs with m = 2, only the 
energy of pseudospin polarized state is independent of 
V2. We therefore systematically control the energy dif- 
ference between Ep=i and Ep^Q by changing V2. The V2- 
depcndence of £'p=i — £'p=o is shown in the inset of Fig. 4 

r , . n j-j^ -j^SLG(n=l) , T^SLG(n=l) 

as a function of ()l/2 = V2 ^ ^2 where V2 

corresponds to the original V2 in the n = 1 LL of SLG. 
This V2-dependence shows systematic change in the size 
dependence of £'p=i — i?p=o at dV2 = —0.002. Since SV2 
is negative at the transition, where Ep^i — i?p=o = hi 
the thermodynamic limit, the pseudospins in the ground 
state at Vn — 2/5 in n ^ 1 LL of SLG is fully polarized in 
large systems. Similar analysis on the size dependence of 
the polarization energy also shows that the pseudospins 
in the ground state of n = 2 LL of BLG are fully polar- 
ized at ^'„ = 2/5. 

The elementally charge excitation energies at Vn ~ 2/b 
are presented in Fig. 5. In the case of rt = LL of SLG, 
the pseudospin (valley) polarized excitation gap Ac is 
calculated supposing pseudospins in the ground state are 
fully polarized, although they are unpolarized in the true 
ground state. The extrapolated value of Ac in the n = 
LL of SLG is then 0.05 e^/(e£p) in a good agreement 
with the previous workj22, In the n ~ \ LL of SLG, the 
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FIG. 6: (Color online) Two-particle correlation functions Qrr' 
of quasihole skyrmion state at Vn ~ 2/5 and quasiparticle 
state at Vn = 2/3 in the n — 1 LL of SLG. r is in units of £b- 
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FIG. 7: The pseudospin (valley) polarized excitation gap Ac 
and the pseudospin (valley) unpolarized excitation gap A3 at 
Vn = 2/3 in n = and 1 LLs of SLG. 



pseudospins in the ground state are fully polarized and 
Ac is 0.059. In BLG, Ac = 0.052, which is also slightly 
larger than Ac in the n = LL of SLG. The largest Ac 
at f„ = 2/5 is obtained in the n = 1 LL of SLG. This 
feature is the same as the case of Vn = 1/3 shown in 
Fig. 3. 

The pseudospin (valley) unpolarized excitation gap A^ 
at Vn ~ 2/5 are also shown in Fig. 5. The A^ from the 
unpolarized ground state in the n = LL of SLG is 0.04 
e'^/{e£B), and the As from the polarized ground state 
in the n 1 LL of SLG 0.035 e^/iels)- These A^ are 
smaller than Ac similarly to the case of Vn = 1/3. In the 
n = 2 LL of BLG, As is larger than Ac. The lowest gap 
in BLG is then given by Ac, which is 0.052 e'^/{eiB)- 

In the n = 1 LL of SLG, the ground state is pseu- 
dospin polarized while the lowest charge excited state 
is pseudospin unpolarized. The pseudospin structure in 
the unpolarized excited states at vi = 2/5 is shown in 
Fig. 2 (b) , which indicates elementally charge excitations 
are characterized by valley skyrmions. The two-particle 
correlation functions gT-T-/(r) ^^i^^i^^ at — 2/5 are pre- 
sented in Fig. 6, which show that the correlation function 
between the electrons in the same valley gKK{f) in the 
unpolarized quasihole state has peak structure around 
the origin while the correlation function between the elec- 
trons in different valleys gKK'if) has the maximum at 
the opposite side on the sphere J^i^ These results indi- 
cate the elementally charge excitations at v„ = 2/5 in the 
n = 1 LL of SLG are valley-pseudospin textures similar 
to skyrmion excitations in the quantum Hall ferromag- 
netic states at Vn = 1 and Vn = 1/3, although their origin 
and detailed properties are not obvious. 

We finally study the valley-polarization in the ground 
state and the elementally excitations at Vn = 2/3. Sim- 
ilarly to the case of the v„ ~ 2/5 fractional quantum 
Hall state, the Vn = 2/3 fractional quantum Hall state 
is mapped onto n — 2 integer quantum Hall state within 



a mean-field analysis of the composite fermion theory. 
Thus the pseudospin unpolarized ground state is ex- 
pected. Indeed, the ground state at ^'„ = 2/3 in the 
n = LL of SLG is pseudospin unpolarized. However, 
in the n = 1 LL of SLG and the n = 2 LL of BLG, 
the size dependence of the polarization energy shows the 
pseudospins are fully polarized in the ground state. 

The charge excitation energies from the ground states 
at i/„ = 2/3 are presented in Fig. 7. The polarized ex- 
citation gap Ac in the n = LL of SLG is calculated 
from the energy difference between the fully polarized 
ground state and the fully polarized charge excited state, 
although the true ground state is the unpolarized state. 
The extrapolated value of Ac in the thermodynamic limit 
is then 0.101 eV(e^B) in the ??, = LL of SLG, which is 
the same as Ac at 1/3 because of the particle-hole 
symmetry of the single component quantum Hall system. 
In the n ~ I LL of SLG, the pseudospins in the ground 
state are polarized and Ac is 0.115 in the thermodynamic 
limit. In BLG, Ac = 0.103, which is also the same as Ac 
at Vn = 1/3. 

The valley unpolarized excitation gap As from the un- 
polarized true ground state in the n = LL of SLG is 
0.08 eV(e^B)- Similarly to the case of 1/ = 1/3, As is 
smaller than Ac. In the n = 1 LL of SLG, we find As 
is 0.10 in the thermodynamic limit, which is also smaller 
than Ac. The valley unpolarized elementally charge ex- 
citations at ^'„ = 2/3 in the n = 1 LL of SLG are not 
skyrmion-like pseudospin-textured state, in contrast to 
the case of = 2/5 in the n = 1 LL of SLG as shown 
in Fig. 6, where short-range correlation gKK is smaller 
than gKK' for = 2/3. 
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TABLE I: The pseudospin (valley) polarization in the ground 
state and the lowest charge excited state, and the unpolarized 
excitation gap As and the polarized excitation gap Ac ex- 
trapolated to the thermodynamic limit at the effective filling 
Vn in the nth LL of single-layer graphene (SLG) and bilayer 
graphene (BLG). 





LL 


ground state excited state 


As 


Ac 


1 SLG 





polarized 


skyrmion 


0.63 




1 SLG 


1 


polarized 


skyrmion 


0.28 




1 SLG 


2 


polarized 


skyrmion 


0.42 




1 BLG 


2 


polarized 


skyrmion 


0.48 




1/3 SLG 





polarized 


skyrmion 


0.05 


0.101 


1/3 SLG 


1 


polarized 


skyrmion 


0.03 


0.115 


1/3 BLG 


2 


polarized 


polarized 




0.103 


2/5 SLG 





unpolarized 


unpolarized 


0.04 


(0.050)^ 


2/5 SLG 


1 


polarized 


skyrmion 


0.035 


0.058 


2/5 BLG 


2 


polarized 


polarized 




0.052 


2/3 SLG 





unpolarized 


unpolarized 


0.08 


(0.101)^ 


2/3 SLG 


1 


polarized 


unpolarized 


0.10 


0.115 


2/3 BLG 


2 


polarized 


polarized 




0.103 



IV. DISCUSSION 

Our DMRG calculation confirms various types of quan- 
tum Hall states in graphene at j/„ = 1, 1/3,2/5 and 2/3 
in the n = and 1 LLs of SLG and in the 7i = 2 LL 
of BLG. These results are summarized in Table I, where 
the pseudospin polarizations for the ground state and the 



lowest charge excited state are listed with the elementally 
charge excitation energies Ac and Ag. The elementally 
charge excitations are obtained by increasing or decreas- 
ing the flux quantum number by 1 for v = 1 and 
1/3,2/5, and by increasing or decreasing the flux quan- 
tum number with adding or removing single electron in 
the system for v = 2/3. We have studied both (a) pseu- 
dospin polarized excitations (Laughlin's quasiholes and 
quasiparticles) and (b) pseudospin unpolarized excita- 
tions (quasihole skyrmions and quasiparticle skyrmions 
at Vn = 1 and 1/3 in the n and 1 LLs of SLG, and 
at Vn = 1 in the n = 2 LL of BLG). 

The activation energies obtained in finite systems are 
extrapolated to the thermodynamic limit, which give the- 
oretical predictions for future experimental studies of the 
fractional quantum Hall states in graphene. Our results 
show that the gaps (Ag ~ O.le^/^s) at 1/3 and 2/3 ef- 
fective fillings in the excited LL in bilayer graphene are 
larger than those in conventional quantum well and sin- 
gle layer graphene. Therefore bilayer graphene is a good 
candidate for future experimental observation of the frac- 
tional quantum Hall effects. 
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